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The quantum states built with the eta paring mechanism i.e., eta pairing states, were first introduced in the 
context of high temperature superconductivity where they were recognized as important example of states al- 
lowing for off-diagonal long-range order (ODLRO). In this paper we describe the structure of the correlations 
present in these states when considered in their momentum representation and we explore the relations between 
the quantum bipartite/multipartite correlations exhibited in k space and the direct lattice superconducting cor- 
relations. In particular, we show how the negativity between paired momentum modes is directly related to 
the ODLRO. Moreover, we investigate the dependence of the block entanglement on the choice of the modes 
forming the block and on the ODLRO; consequently we determine the multipartite content of the entangle- 
ment through the evaluation of the generalized Meyer Wallach measure in the direct and reciprocal lattice. The 
determination of the persistency of entanglement shows how the network of correlations depicted exhibits a 
self-similar structure which is robust with respect to "local" measurements. Finally, we recognize how a rela- 
tion between the momentum- space quantum correlations and the ODLRO can be established even in the case of 
BCS states. 
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I. INTRODUCTION 

The eta pairing mechanism for electrons was first intro- 
duced by Yang 1 1] in the context of high- temperature super- 
conductivity. Generally speaking, the relevance of the states 
built by means of the eta pairing mechanism, i.e. eta pairing 
states, comes from the fact that they are exact eigenstates - in 
some cases actual ground- states - of several relevant models 
in different fields of condensed matter physics. 
A first example of a model having the eta pairing states as 
eigenstates is the Hubbard one. It can be proven 1 1] that in this 
case these eigenstates cannot be the ground- state since one 
can build states with different symmetry having a lower en- 
ergy. Nevertheless, the eta pairing mechanism allows for the 
formation of Cooper pairs located on a site, rather than sepa- 
rated by a finite coherence length, as in the case of the standard 
Bardeen, Cooper Schrieffer scheme for superconductivity | 2]. 
In this context these eigenstates play a prominent role since 
they display off-diagonal long-range order (ODLRO) 1 3]. The 
latter is a peculiar kind of nonlocal correlation that survive 
in the thermodynamical limit. It has been shown |4] that 
ODLRO implies Meissner effect and flux quantization, which 
are both distinctive features of superconducting systems. 

On the other hand, for specific regimes of parameters, the 
eta pairing states turn out to be the ground- state of different 
extensions of the Hubbard model |[5[_0, S Isl, Hi- In partic- 
ular, in the bond-charge extension |9] a state characterized 
by the eta pairing mechanism belongs to the lowest-energy 
sector even at finite positive values of the Coulomb interac- 
tion. Also, by mapping the eta pairing state into the spin-^ 
language |8], one can recognize it as the ground state of the 
isotropic Heisenberg ferromagnet, as well as of the Heisen- 
berg antiferromagnet when the anisotropy parameter A = — 1 
LlO]. 

Due to their relevance, the eta pairing states are natural can- 
didates for the applications of the tools and the paradigms de- 



veloped within the framework of quantum information theory 
(QIT) [11, 12]. In this context, a central concept is that of 
entanglement i.e., the quantum correlations among or within 
subsystems constituting a system. A significant amount of re- 
sults in QIT aims at the classification and quantification of the 
existing correlations, quantum or classical, in a given quan- 
tum state. The Von Neumann entropy, for instance, measures 
the quantum correlations of a subsystem with the rest of the 
system. The quantum mutual information is a measure of the 
total (quantum and classical) correlations between two sub- 
systems, while the negativity quantifies just the quantum cor- 
relations between the latter. One can thus gain a deep insight 
about the structure of the existing network of correlations in 
a complex system by evaluating and comparing appropriate 
entanglement (correlation) measures once the sets of different 
partitions of the system into subsystems are given. 

In the direct lattice picture, Zanardi and Wang ifTsIl were 
the first to apply the above scheme to the eta pairing states 
by analyzing the entanglement between two sites and show- 
ing that the latter vanishes in the TDL. The entanglement be- 
tween a block of sites and the rest of the lattice was analyzed 
in lUn, H \l5] where it was shown that it scales logarithmi- 
cally with the size of the block and how it is connected with 
the ODLRO. A through analysis of the entanglement in the 
direct lattice picture for pure states and mixtures of eta pair- 
ing states was also carried on in lfl6ll . where it was pointed 
out that while the two- site entanglement vanishes in the ther- 
modynamical limit (TDL), the two-site classical correlations 
are still present; when measured by the quantum mutual in- 
formation, the latter were recognized in LIT ] to coincide with 
the ODLRO. Finally, in O H it was investigated how the 
quantum (total) correlations of ground- states characterized by 
the presence of the eta paring mechanism behave at quantum 
phase transitions. The analysis allows one to classify the latter 
in terms of two- site or multipartite quantum correlations. 

In this paper we aim at studying the entanglement proper- 
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ties of the eta pairing states in their momentum representation. 
This study on one hand provides a complementary and richer 
picture of the underlying correlations structure; on the other 
hand it allows one to explore the relations between the quan- 
tum bipartite and multipartite correlations exhibited in k space 
and the superconducting correlations typical of the direct lat- 
tice picture. In particular, we evaluate appropriate measures 
of correlations considering different choices of the elementary 
subsystems (single, paired and unpaired modes) and we find 
how the negativity between paired modes is directly related 
to the ODLRO present in these states. Furthermore, the exact 
diagonalization of the reduced density matrix associated with 
an arbitrary set of momenta will allow us to study the scaling 
of the block entropy with the size of the block and to see how 
it depends on the selection of the modes forming the block 
and on the ODLRO. These results will also make it possible 
to to compare the multipartite entanglement content in the di- 
rect and reciprocal lattice picture through the evaluation of the 
Meyer Wallach measure [19, 20, 21]. The determination of the 
persistency of entanglement ll22ll will show how the network 
of correlations depicted exhibits a self- similar structure which 
is robust with respect to "local" measurements, i.e. measure- 
ments of single or paired momentum modes [23] . 

In order to compare the results obtained in the eta paring 
case we finally study the BCS states. Some entanglement 
properties of these states were studied in [2^, [25l [3 l27h . 
Here, by resorting to the Green's function language developed 
in 1I27I1 we show that while these states exhibit a simpler struc- 
ture of correlations in momentum space this structure can be 
directly linked to their ODLRO. 

The paper is organized as follows: we start by introduc- 
ing the /c-space representation of the eta pairing state then, in 
section HTbI we analyze the correlations between single modes 
and pairs of modes. Section [ncl is devoted to the study of the 
block entropy, while in Section HTPl we review the Meyer Wal- 
lach measure and we use it to investigate the multipartite con- 
tent of correlations in the eta pairing states. We conclude the 
analysis of the correlation properties of the latter in Section 
III El where we analyze the persistency of the entanglement. 
Finally, in Section [nil we investigate the connection between 
entanglement in momentum space and ODLRO in BCS states. 

II. THE ETA PAIRING STATE 

A. Preliminaries 

The eta paring states are built through the action of the so- 
called eta operator that, in direct lattice picture, is written as 

^^ = E4T4i = E^r (1) 

Here rjj = cj^cj^ is written in terms of cj^ and cj- i.e., the 
fermionic creation operators on the site / of a one-dimensional 
chain of length L; a G {j, |} is the spin label and a de- 
notes its opposite. When acting on the vacuum of the lat- 
tice it creates a pair of particles fully delocalized 



over a chain of sites of length L. By using the Fourier trans- 
form a|. = YliLi ^^^^^^la ^^^^ 4a obtains the 
/c-space representation of the eta operator: 

L-l L-1 

j=0 j=0 

where 7}^ = a^^^a^ Each kj = 27rj/L now labels 
one of L momentum modes, whose local basis is Bkj = 
{|0)/c^.,| T)/c,-,| i)kjA U)kj}- When acting on the vacuum 
|vac)x = ®f=o the eta operator creates a pair of parti- 
cles fully delocalized over the whole momentum space. The 
way in which the delocalization is performed in the two repre- 
sentations is fundamentally different. In the direct lattice each 
of the r]j operators acts creating a pair of particles (T, j) local- 
ized on the site In the k representation, the delocalization is 
performed already at the level of each fjl ; the latter acts cre- 
ating pairs of particles of the type (cr-kj , o^kj ), thus involving 
the pair of modes {—kj, kj ) . 

The generic eta pairing state is obtained by the creation of 
Nd pairs of particles: 



\nL,N,)) = ^i^^(r,t)^<<|vac) . (3) 

We note that the states considered here are a particular case 
of the general family of eta pairing states that can be built by 
operators of the form r/J^ = Ylf=i ^^^^^^^i- the /c-space 
picture these operators create pairs with momentum equal to 
(j). One can see that the structure of correlations of the states 
generated by the action of r^J^ does not depend on the partic- 
ular choice of the value of ^; in the following we will thus 
choose ^ = 0. 

Our analysis is based on the evaluation of the (quantum) cor- 
relations between different subsystems i.e., set of momentum 
modes. To this end we consider the following measures of cor- 
relations. The Von Neumann entropy will be used to measure 
the quantum correlations between a subset A of momentum 
modes and the rest of the system when the latter is in a pure 
state. Its definition is based on the reduced density matrix pA- 

S{pa) = tr {pA log2 pa) = - ^ Ai log2 (4) 

i 

where {A^} are the eigenvalues of pA- In order to measure 
the total i.e, quantum and classical, correlations between two 
generic subsystems A and B (sets of momentum modes), we 
use the quantum mutual information lO, [28, 29]. This mea- 
sure is defined in terms of the system's and subsystems' den- 
sity matrices pab, Pa and pb' 

Iab = S{pa) + 5(pb) - S{pab)' (5) 

The quantum correlations between two generic subsystems A 
and B will be quantified by the negativity 1 30] 

A/'(Pab) = (||pI1j||i-1)/2; (6) 
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where p^j^ is the partial transposition with respect to the sub- 
system A appHed on pab, and ||0||i = TrVOW is the 
trace norm of the operator O. All the above functionals prop- 
erly capture the bipartite correlations between two subsys- 
tems. For the multipartite correlations we will make use of 
the Meyer Wallach measure see section HTPl for 

the definition. 

As mentioned in the introduction, the eta pairing states were 
first recognized as relevant since they allow for superconduct- 
ing correlations i.e., ODLRO, which for these states are de- 
fined as: 



lim {r]^.r]j 



i+r] 



nd(l - rid) • 



(7) 



where j,j-\-r label two sites of the lattice at distance r. A 
first bridge between ODLRO and a measures of correlations 
defined in the context of quantum information theory was es- 
tablished in 1 17] where the analysis carried on in the direct 
lattice framework allowed to recognize that these correlations 
are simply proportional to Tjj-\-r i.e., the mutual information 
between two generic sites j^j-\-r, where r is an arbitrary dis- 
tance . 



B. Correlations in /c-space: single and pairs of modes. 

In this section we begin to analyze the structure of corre- 
lations of the eta pairing states at the local level. The calcu- 
lations involved in the study of both the single mode and the 
pairs of modes will lead in the following sections to the gen- 
eralization to the case of block of modes. 



1. Single mode 

We start the analysis by evaluating the reduced density ma- 
trix of the single generic mode kj . The calculations are 
instructive since some of the features that hold for the single 
mode case will be useful for the determination of the reduced 
density matrix of the set of modes we consider in the follow- 
ing sections. 

It turns out that the only non- vanishing matrix elements of pk^ 
are the diagonal ones. The latter are the expectation values of 
the following projectors: 

«fe,i«/c,T4,T4,i = |0)/c,/c,(0|, 
4,T^fc.T«fc,i4,i = I T)/c,/c,(T I, 

^kji^kjiCLkj^CLl,^ = I i)/c,-/c,-(i I, 

^k,^^k,i^k,iak,^ = I M)k,k,{[] I . 

Indeed, the state \^{L^Nd)) is built through the action of 
superpositions of operators f}\, : every time the mode kj is 
occupied, the —kj one is changed correspondingly. There- 
fore, since the off-diagonal elements correspond to projec- 
tors that represent processes that, acting locally onto the kj 



mode, leave unchanged the —kj one they have zero expec- 
tation value. Accordingly, the diagonal single-mode reduced 
density matrix reads: 



Pk^ =6:mg{AA',AB,AB,BB'} 



(8) 



where 



A' = ^] ^ 
B = ^ I L^oo 



B' = 



a = 1 — Ud 



b = l-a = nd. (9) 



L-l 



Capital letters refer to finite-size expressions while lower 
cases refer to their asymptotic expression for L oo. We 
can now evaluate as first measure of correlations the Von Neu- 
mann entropy of pk^ : 



Sk, = -2 (a log a + 6 log b) 



(10) 



which gives the amount of total (quantum) correlations that 
the single generic mode has with the rest of the system. Just 
as much as in the direct lattice picture (|[i3,[ill]), the correla- 
tions are directly determined by the filling Ud = n/2. and they 
reach their maximal value in correspondence with the half fill- 
ing case n = 1. 



2. Two modes 

We now consider the correlations of the subsystem consti- 
tuted by two generic modes {ki^kj). We have to distinguish 
between two cases. For ki ^ —kj the reduced density-matrix 
^ is diagonal with respect to the local basis S/c . ®Bkj ; in- 
deed, as described for the single mode case, the off-diagonal 
elements correspond to expectation values of projectors that 
change the state of the modes ki, kj and do not affect the 
modes —ki,—kj \ they hence have zero expectation value. For 
large L, the eigenvalues are a'^""^", a = 0, . . . , 4, each one 
appearing with multiplicity given by rrioc — (^) . This scheme 
can be generalized to a higher number modes, as we shall see 
in section HTCl The case ki = —kj has to be treated sepa- 
rately. The support of the reduced density-matrix is the sub- 
space spanned by 

= {|o,o),-, I ra),, I i,T)„ I n,n>,} , (id 

where |Q^,/3)j = \oi) -■k-\tS)]^.. Indeed, the sole states that 
can be built by the action of the r)^ operators belong to this 
subspace. Thus the reduced density matrix of the subsystem 
(— /cj, kj) has just a 4 x 4 nonzero sub-block relative to the 
subspace spanned by (fTTI) : 



/a^ 

ab ab 

ab ab 

V 6^ 



(12) 
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whose diagonal form is diag{a^, 2a6, 0, 6^}. 

The density-matrices just determined allow us to calculate 
i) the quantum correlations between the pair of modes with 
the rest of the system ii) the total and quantum correlations 
between any pair of modes (/c^, kj). We have that: 



—2 (a log a + 6 log 6 - 
—4 (a log a + 6 log 6) 

Ski + 2a6 ki = —kj 

^ ' 

ab/3 ki — —kj 
ki kj 



ah) ki = 



kj 



kj 



(13) 



(14) 



(15) 



where: Sk^^kj^ the Von Neumann entropy of ^ , mea- 
sures the quantum correlations between the pair of modes and 
the rest of the systems. Tk^^kj, the mutual information, and 
Mki.kj , the negativity, measure the total and quantum correla- 
tions between the modes. 



3. Paired modes 

In order to understand the /c-space network of correlations, 
at the local (and multipartite) level, we turn to the natural 
extension of the above calculations by considering the four 
modes {—kj^kj)^{—ki^ki). This will allow us, on one hand, 
to identify as elementary subsystem the generic pair of modes 
{—ki^ki) and to study the nature of the correlations between 
such kind of subsystems, i.e. two-pair correlations. On the 
other hand, the strategy at the basis of the calculations can be 
generalized to an arbitrary number of pairs of modes, thus en- 
abling the evaluation of the entropy of a block of modes with 
respect to the rest of the system (block entropy). 

We first must determine the reduced density-matrix 
fc -fc fc ) • turns out that the latter has support only on 
thesubspacespannedby _^ = ^-kiM'^^-kj,kj, 
see (fTTt . In fact, the only 4-modes projectors that have 
a non vanishing expectation value are those that preserve the 
number a of {cr-kj , d^kj ) pairs of particles (as in the two-mode 
case). Moreover, the expectation values are all equal for a 
given a and can be straightforwardly determined as functions 
of a. For large L we have: 



a i 
ah 



= aH 



ah'' 



(16) 



i.e., (P^) = a^~"6" is the expectation value of the pair- 
preserving projector between local states containing a pairs. 
The fact that the expectation values of the elements involv- 
ing a given number of pairs are all equal allows for a simple 
expression for the of diagonal form of the non vanishing sub- 
block of p. 



a fixed number of pairs. The latter is, in fact, equal to the 
number of ways one can set a pairs in 4 modes. 

We have now all the ingredients to derive the following cor- 
relations measures: 



5, 



( ki,ki),{ kj,kj) 



-4 {a log a + 6 log 6 + 1 /2a(l - a) • 
• [4 + a(l-a)(31og2 3-5)]} (18) 
2a(l - a) [2 + a(l - a)(3 log2 3 - 5)] 

(19) 

(20) 



Furthermore, the above evaluated reduced density matrices al- 
low us to determine that the generic single mode kj is un- 
corrected with respect to any pair of modes {—ki^ki), i.e. 

^kj,{-ki,ki) = 0- 




diagjmoa^, mia^6, m2a^6^, rrisah^^ m4h^} 

(17) 

The coefficients rua = (^) are given by the dimension of 
the square sub-matrix of p(^_ki,ki-kj,kj) that corresponds to 



FIG. 1: Left: schematic representation of the correlations of the 
generic single mode —ki and of the paired modes {—ki,ki). QC 
stands for quantum correlations, CC for classical correlations, while 
QMP stands for quantum multipartite. Right: two opposite ways of 
grouping the modes in blocks i.e., paired (Adi) and unpaired (AD2) 
modes. 



4. Discussion 

The above /c-space calculations suggest the following sce- 
nario sketched in figure[T] The quantum correlation of the sin- 
gle mode kj with the rest of the system have both a two-mode 
and a multipartite nature. The mode has two-point quantum 
correlations only with the mode {—kj), see ([T5]) , while all the 
other two-mode (total) correlations are zero: when ki 7^ —kj, 

^kj,ki = 0. 

The fact that S-kj,kj 7^ 0, see ([13]), implies the existence 
of a quantum multipartite contribute to the single-mode corre- 
lations. The nature of such a contribute is better characterized 
by the fact that the single mode is uncorrected with any pair 
• of modes ki^ki 7^ —kj, i.e., Tk^^kiki = 0- This means that we 
are in presence of n-way entanglement, with n > 3. All the 
considered measures reach their maximum at a = 1/2; the 
latter condition is satisfied when Ud = n/2 = 1/2. We now 
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focus our attention on the elementary subsystems constituted 
by the pairs of modes (— /cj, kj). These are quantum corre- 
lated with the rest of the system ([T3]) and in order to describe 
the multipartite correlations related to them we can introduce 
the notion of multi-pair entanglement; indeed, the mutual in- 
formation ([T9l) indicates the existence of two-pair correlations 
which are of classical nature only, since at the same time, the 
two-pair quantum correlations, as measured by the negativ- 
ity (|2Q1) , are zero. This kind of distribution of multipartite 
correlations in terms of two-pair classical ones is the same 
encountered in the direct lattice picture where a the level of 
two-sites the correlations are just classical and proportional to 
the ODLRO. 

It is now important to explore the nature of some of the 
described correlations and highlight their relations to mean- 
ingful physical quantities. In particular it turns out that it is 
possible to establish a direct connection of the quantum cor- 
relations between two paired modes and the ODLRO. Indeed, 
the negativity M-kj^kj is calculated on the basis of the off - 
diagonal elements of the non- vanishing part of the reduced 
density matrix of the paired modes ([T2l) . These elements co- 
incide with expectation values {f}l,f}-kj) = (^1/^ fjkj ) and by 
writing them in terms of direct lattice operators one has: 

A4,,-fe, = {4,'?-fe,)/3 = ^ Y^^riUra) + (21) 

The sum has L{L — 1) non vanishing terms and we see that 
the quantum correlations between paired momentum modes 
are given by the average of the pairing correlations in the di- 
rect lattice. Since the latter do not depend on the lattice sites' 
indexes m and are all equal to ^ we the above equation 
establishes a direct connection between the quantum corre- 
lations between paired momentum modes and the sufficient 
condition for superconductivity i.e., ODLRO. As a final re- 
mark we note that the ODLRO (|7]) explicitly appears in the 
expression of other correlation measures, in particular those 
indicating, as previously discussed, the presence of multi-pair 
entanglement ([T3l) and its reflection at the level of two-pairs 
correlations([T9b. In the following sections we will see how 
this multi-pair nature of quantum correlations determines the 
behaviour of the block entanglement and of the related mea- 
sure of multipartite entanglement and how these can be again 
expressed in terms of ODLRO. 

C. Block entropy 

We now deepen our study of the structure of the correla- 
tions in the eta-paring state {^{Nd)) by considering how the 
subsystems constituted by blocks of modes are entangled with 
the rest of the system: block entropy. 

In the direct lattice picture the eta-pairing state can be mapped, 
through a particle-hole transformation |8], onto the ground 
state of the isotropic ferromagnetic Heisenberg model. In this 
framework, by taking advantage of the permutational invari- 
ance of such states, one can evaluate the entropy of blocks of 
sites 



Here, in the /c-space representation, the structure of correla- 
tions between blocks of modes is more complex and allows 
for a richer picture: as we shall see in the next section, the 
general form of the block entropy depends on the particular 
choice of the set of modes constituting the block. 

1. Unpaired modes 

In this section we evaluate the block entropy between a set 
Ad^ of Di single modes and the rest of the system such that 
Mkj ^ Ad^ ^ —kj^Aoi- this choice of set of modes 
the reduced density matrix p^i of the block is diagonal since 
the off-diagonal elements correspond to projectors that act 
only on the /cj-mode while leaving the —kj mode unchanged 
and therefore have zero expectation value. The generic diago- 
nal element is the expectation value of the product of I^i local 
projectors Pfc^. e {|0)(0|, | T)(T 1,1 1,1 Ti)(Ti I}- The 
result can be written in terms of the number of modes and the 
total number of fermions M involved by the projectors: 

(L-2Di\ L ^ OO 

{K')= /l^ a^^^-n^. (22) 

There are {^^^) projectors that involve M particles and they 
all have the same above expectation value. Thus the von Neu- 
mann entropy of the block corresponding to the set Ad i can 
be written as 

= -2Di (a log a + Mog6) = DiSk^ . (23) 

Indeed, this result reflects the factorization of the density ma- 
trix of the block 

pDi = ^kjeAo^Pkj , (24) 

thus extending the result ([T3l) for the single mode kj . Equation 
(l24b shows how the modes belonging to Adi are completely 
uncorrected among each other. The correlations of the block 
are given by i) multipartite correlations between each mode 
and the rest of the system ii) the two-mode correlations: the 
boundary of the block cuts the links between each kj G Adi 
and its complementary —kj ^ Adi,^^^ figure [T] 

2. Paired modes 

In this section we consider a different situation which is 
opposite, in some sense, with respect to the previous one. We 
evaluate the block entropy between a set Ad2 of D2 single 
modes and the rest of the system such that \/kj ^ Ad2 ^ 
-kj e Ad2- 

The results obtained in section IIIBI for one and two pairs 
of modes can be extended to an arbitrary number D2/2 of 
pairs of modes {—kj^ kj). In particular, the expressions ([T6b 
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obtained for the two pairs of modes can be easily general- 
ized. The projectors t^^t have nonvanishing expectation 
value are those that conserve the number of pairs of particles 
{(T-kj , ^kj ) and that represent, as explained above, processes 
that "coherently" affect the (— /cj, kj) pair. 
The density matrix relative to 1^2 /2 pairs of modes is block- 
diagonal. Each block can be labeled by the fixed number of 
pairs a G [0,1^2] of (cr_/e^ , a/e^ ) particles involved by 
The elements of a given block are all equal and their explicit 
expression is 



(L-D2\ 
\Nd-a) 



oo 



(25) 



The dimension of the block is rria = (^^). i-^-. it corre- 
sponds to the number of ways one can place a pairs of parti- 
cles {(T-kj , ^kj ) in D2/2 pairs of modes. Each block labeled 
by a has rria — 1 null eigenvalues. In the large-L limit the non- 
vanishing part of the spectrum of the reduced density matrix 
is then given by the following set of D2 eigenvalues: 

{moa^% mia^2"^6, . . . , mo^-iab^^'^ ^rno^b^^} • (26) 

Since b = 1 — a we see that the eigenvalues follow a binomial 
distribution and the Von Neumann entropy thus reads 



D2 

E 

M=0 



D2 

M 



D2 

M 



(27) 



For large values of D2 the latter expression was evaluated 
in 1321 and it has the following asymptotic expression: 



1 



1 



(28) 



that holds for a6 > 0. 



3. Paired and un paired modes 

We have seen that the choice of the set of modes crucially 
determines the (asymptotic) behavior of the block entropy. In 
particular, the block entropy of the set of modes Adi , which 
does not contain any pair (— /cj, kj), grows linearly with the 
number of modes ([23l) . On the contrary, in the opposite situ- 
ation, the block entropy of the set of modes grows only 
logarithmically. 

We note that results (|27]) and (l28l) are only formally equiva- 
lent to the one obtained in the spin-models context 1 15] since 
the latter refer to the direct lattice picture and they take into 
account completely different subsystems (blocks of sites). In- 
deed, while in the direct lattice the choice of the block of 
modes is unambiguous here, in the /c-space, the block entropy 
behaves logarithmically only for a specific choice of the set 
of modes, i.e, Ad2' Furthermore, in order to give a complete 
description of the block entropy one has to consider the inter- 
mediate picture of a block composed hy D = Di^D2 modes 



belonging to the ^Qi Adi+D2 = ^Di\JAd2^ where Di is the 
number of unpaired modes kj G Adi and 1^2/2 is the number 
of the paired modes {—kj^kj) G Ad2- The generic nonzero 
element of the reduced density matrix of the block is given by 
the expectation value of projectors of the type 

¥d{Di, F, a) = ¥d, {F - 2a) Pi^, {a) . (29) 

The latter are characterized by a fixed number of fermions 
F = M -\- 2a: M is the number of single particles involved 
by the Di single mode projectors P/e^. G {|0)(0|, | T)(T |, | i 
)(i U Ti)(Ti |}' while a is the number of pairs of particles 
{(T-kj , o-kj ) involved by the D2 projectors relative to the pairs 
of modes. 

It turns out that the generic expectation value can be written as 
the product of expectation values ofFDi{F—2a) and Pi:)2 (a), 
i.e. 

{¥D{Di,F,a)) = {¥DAF-2a))-{FD-DM)) = 



I) 



L-2D 



-D + Di 



(30) 

This expression can be used to recognize that, for large L the 
contributions of paired and unpaired modes can be factorized; 
indeed, for finite L, the reduced density matrix of the block 
can then be written as: 



PD = pDi ^ PD2 



(31) 



and this implies that for any set Adi the unpaired modes are 
uncorrected with respect to any finite number D2 of paired 
modes {—kj^ kj). If we now let L take arbitrarily large values 
formulas ([3Ub can be simply expressed in terms of a, 6 = 1 — a 
as 



^D{Di,F,a)) = a^+^^-^+"6 



(a 



D — Di—a-La.\ 



(32) 



and one can use this relation to recognize again that [see 
([T4b l for L large enough the unpaired modes are uncorrelated 
among each other, i.e. 



PD^ 



Pki 



(33) 



Since for fixed Di and D2 the structure of p^^ and pD2 is 
the same analyzed in the previous sections the evaluation of 
the eigenvalues pD follows along the same lines and the Von 
Neumann entropy of the block of modes in the set Adi+D2 
can thus be simply expressed as the sum of the block entropies 
pertaining to Ad^ and Ad2 ' 



Sd^+D2 



-2Di (a log a + 6 log 6) 



D2 

E 

M=0 



M 



'6^ log 



D2 
M 



D2-MrM 



(34) 
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When D2 is sufficiently large on can again approximate the 
binomial distribution and obtain: 

^1^1+1^2 - -2Di{a\oga^b\ogb)^^\ogD2 . (35) 



one can further generalize and extend the measure by consid- 
ering principal subsystems Am composed by m elementary 
subsystems S. Here the average must be extended to all the 
= n\/m\{n — m)\ possible inequivalent choices of Am 
and it can be written as 



4. Discussion 

The above results suggest the following interpretation in 
terms of two-point and multipartite entanglement. On one 
hand the linear growth of the block entropy Sdi has two main 
contributions: i) the two-mode quantum correlations of each 
mode kj e Adi with the corresponding mode —kj ^ Adi', 
a) the multipartite correlations that a single mode shares with 
the rest of the system. On the other hand the logarithmic 
growth cannot be described in terms of two-pair quantum 
correlations: since M(^-kj^kj),{-ki,ki) = 0^ the growth of 
reflects the multipartite (multi-pair) nature of the correlations 
of the subsystem Ad2 with the rest of the system. In the inter- 
mediate case Adi U Ad2^ these two effects sum and Sdi+D2 
accounts for both the two-pair and the multi-pair contributes 
to correlations. We finally observe that the binomial distribu- 
tion (|26b , that corresponds to the eigenvalues of the density 
matrix relative to block of pair of modes, contains the link 
with the ODLRO. Indeed, the variance of the distribution is 
just D2nd{l — rid), thus to a higher amount of ODRLO it cor- 
responds a broader probability distribution and consequently 
a higher value of the block entropy (l27l) . The above discussed 
multi-pair contribution to the entanglement present in the mo- 
mentum picture is thus directly related to the superconducting 
correlations. An analogue relation was find in the direct lat- 
tice picture m id H [S m where the entropy of block of 
sites was considered and an analogous binomial distribution 
for the eigenvalues of the reduced density matrix was found. 
In the following we will compare the two pictures on the basis 
of the Q-measure of entanglement. 



D. Q-measure of entanglement 

We use the above results to compute the Meyer- Wallach 
measure of multipartite entanglement Qdip)). The latter was 
first introduced in [19] for multi-qubit pure states. In [20] it 
was shown how Q can be simply expressed as the average 
linear entropy of the single qubits: 



Q{m = 2{ 1 



k=l / 



(36) 



here pk is the reduced density matrix of the k-th of n qubits. 
This expression of Q allowed in (2?] for two different kinds 
of generalization. On one hand, one can extend the measure 
to the case of multi-qudit states, i.e. to quantum systems that 
are composed by n identical elementary subsystems S. Each 
of the latter lives in Hd (dim(7Yd) = d), while the whole 
system lives in (C^)^"". In this case, Q measures the aver- 
age linear entropy of the subsystems S. On the other hand. 



1 



1 



(37) 



where each i = {ii, . . . , z^} is a subset of indexes identify- 
ing the m elementary subsystems S composing a given Am, 
the latter being characterized by its reduced density matrix p^. 
The factor d^ / {d^ — 1) normalizes the measure to one. 

Qm,d is thus a multipartite entanglement measure based on 
the average of a bipartite one (linear entropy): it quantifies the 
average entanglement between blocks of qudits and the rest 
of the system. Anyway, only as the size m of the blocks in- 
creases, Qm,d becomes really sensitive to correlations of an 
increasing global nature |2Uf33l]. This property is in general 
difficult to exploit since it requires the evaluation of the re- 
duced density matrices of blocks of arbitrary dimension. As 
we have seen in the previous sections this turns out to be pos- 
sible in /c-space. 



1. Q-measure in k-space 

In this section we analyze the multipartite correlations in k- 
space as they can be described by (5m,d- We have to specialize 
([37I) to the /c-space case, d = 4 is the dimension of the single- 
mode Hilbert space. Each principal subsystem Ad=Di-\-D2 is 
composed by an even number D of modes; Di are unpaired 
modes ki while D2/2 are pairs of modes {—kj^kj). For a 
given L, one has (^) inequivalent choices. In order to eval- 
uate the density matrix of the generic Sd one has to take 
into account its composition in terms of paired and unpaired 
modes. As we have seen in section ITTCl in the finite size case, 
Pd can be written as a tensor product PD2 ^ PDi but the un- 
paired modes are not uncorrected, i.e. PDi 02 Pki- For 
fixed D2 the various po have the same spectrum and since 
Tv{p'jj) = Tr(p|)2 )Tr(p|) J, Q can eventually be written as 



Qd,4 



where 



f{D2 



4^-1 



D 



-1 D 



J2 fiD2)npijTvipij 

(38) 



D2=0 



U^-\l - 2z) U^!,'-\L - 2D, - 2j) 



(39) 



1^1! D2\ 

is the number of equivalent partitions of D into Di = D — D2 
single modes and D2/2 pairs of modes; 

M 



21^1 /, 



a=0 



X (L-D2\-i^ 
\ \Nd-a) 



a 



(40) 
(41) 
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Formula ([38]) can be numerically evaluated for L large; in fig- 
ure [21 we plot (5d,4 for different values of I) as a function of 
Nd/L. As D grows, the measure, which is normalized to one, 
rapidly saturates to its maximal value for any filling. This re- 
sult confirms the analysis of the structure of correlations car- 
ried on in the previous sections. Indeed, in the /c-space pic- 
ture, multipartite entanglement is the dominant feature of the 
eta pairing state. 




0.2 0.4 0.6 0.8 1 ^ 



FIG. 2: (Color online) Plot of the measure of multipartite entan- 
glement Qda in /c-space for different sizes of the elementary sub- 
system; the size varies form bottom to top: D = 1,2,4,8,16 
(L = 1000) 



2. Q-measure in direct lattice 

We conclude by considering the extension of the above 
analysis to the direct lattice picture. As already mentioned, 
through a particle-hole transformation it is possible to trans- 
late the eta-pairing states into the spin language. At variance 
with the /c-space picture, where the idea of considering the 
pairs of modes (— /cj, kj) is natural and necessary, in the di- 
rect lattice the choice of the set of elementary subsystems (the 
sites of the chain) constituting the generic principal subsystem 
(the block of sites) is unambiguous. The local Hilbert space 
has effective dimension d = 2, since each site can be either 
empty or doubly occupied. Due to permutational symmetry, 
each block of sites has the same po and, consequently, the av- 
erage linear entropy coincides with the linear entropy of one 
block. In figure[3]it is shown how, as D grows, Qd,2 asymp- 
totically approaches its maximal value. This result confirms 
that in the direct lattice picture the multipartite entanglement 
does play a fundamental role in the eta-pairing states 1 16i.ll7il. 

It is now interesting to compare the structure of correlations 
scenario emerging from the two frameworks: direct and recip- 
rocal lattice. 

In the direct lattice there are no two-point quantum correla- 
tions: the two- site concurrence is zero when L ^ oo for 
any pair of sites, thus implying a vanishing entanglement ra- 
tio nl\\ . This is a signature that the nature of the correlations 
in this picture is basically multipartite. Furthermore, the two- 
point correlations, i.e. the mutual information between two 



Qd,2 



1 100 




0.2 0.4 0.6 0.8 1 L 



FIG. 3: (Color online) Plot of the measure of multipartite entan- 
glement Qd,2 in direct lattice picture for different sizes of the el- 
ementary subsystem; the size varies form bottom to top: D = 
4,8,16,32,100 (L = 1000) 



sites, since they do not depend on the distance between the 
sites are uniformly spread over all the chain; as previously 
mentioned their value is proportional to the ODLRO. 

In the /c-space picture the only "two-point" quantum corre- 
lations -as measured by negativity- are those shared by each 
mode kj with its "complementary" —kj and we have seen how 
they are directly connected with ODLRO. Furthermore, the 
total two-point correlations are localized between these two 
complementary modes (X/e-^/c^ 7^ ki = —kj). As 

discussed in section [irB 41 it is possible even from this coarse- 
grained level of description, infer the relevance of the multi- 
partite kind of entanglement. 

We can thus pass to analyze the main differences at the level of 
multipartite correlations, as measured by Q. The comparison 
of figure [2] and figure [3] points out that in general, for a given 
size D of the principal subsystem (i.e. block of sites/modes), 
the value of Q is higher in the /c-space case. Indeed, in the 
direct lattice for fixed D, Qd,2 is proportional to the linear 
entropy 

In the /c-space, Qd,4^^ the average of terms that take into ac- 
count the composition of the blocks in paired and unpaired 
modes and it is of course greater than the least of these terms. 
The latter corresponds to the block composed by 1^/2 paired 
modes and formally coincides with ([42l) , see (I39II41I) . These 
considerations suggest that while the dominant feature of the 
eta pairing state is the multipartite entanglement, this kind of 
correlations play a major role in k- space representation. The 
difference between the two representations manifests itself in 
the distinct asymptotic behavior of the block entropy in the 
two pictures: in direct space Sd can grow at most logarithmi- 
cally with the size of the block, while in /c-space it can grow 
linearly ([35]). 
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E. Persistency of entanglement 

We conclude our analysis of the structure of correlations in 
the eta pairing states by briefly discussing -both in direct and 
in reciprocal lattice- the operational effort required to destroy 
all the entanglement present in the system. To this end com- 
pute the persistency of entanglement Ve that was introduced 
in [22] to test the strength of quantum correlations present in 
a state in view of its use for quantum information protocols. 
For a L-qudit pure state Ve is defined as the minimum num- 
ber of local measurements that reduces the entangled state to 
a product state of the L qudits (i.e. that disentangles the state). 

In the direct lattice we have seen that the local state is a 
qubit, since the j-th site can only be either empty or doubly 
occupied. Due to the structure of the eta pairing states, the 
effect of a single local measure Pq, = ^i^jti (g) \a)jj{a\, 
a e {0, 11} is to reduce the number of sites over which the 
eta pairs can be delocalized from L to L — 1 and to generate a 
new eta pairing state with either A^c^ — 1 or Nd pairs, i.e., 

¥o\^{L,Nd))/po = |O),,(O|0|^(L-1,7V,)) (43) 

(44) 

where = Tr{¥a\'^ {L, Nd)){'^{L, Nd)\). The eta pairing 
states thus display a self- similar behavior under local mea- 
surements. Moreover, they turn out to be robust to noise since 
for Nd < L/2 [Nd > L/2] one needs to perform at least 
Nd [L — Nd] measures in order to destroy all the quantum 
correlations present in the state, i.e.. 

Indeed, it easily seen that the minimum number of measures 
needed to factorize {"^{Nd)) is reached by repeating Nd mea- 
sures of type (|43]) , if A^^i < L/2, or L — Nd measures of type 
ei, if Nd > L/2. 

Despite the fact that in momentum space, being the local 
space four-fold, the number of possible local measures is four 
Fa = <S)i^jti \a)jj{a\, a e {0,t,i,Ti}, the situation 
is similar to the direct lattice one, since the more convenient 
sequence of measures to factorize the state consists again in 
performing just measures of type when Nd < L/2 or of 
type Po when Nd < L/2. In /c-space, however, the effect of 
one of these local measures is to reduce the available space 
for the delocalization of the eta pairs from L to L — 2 and to 
generate an eta pairing state, consisting of either Nd or Nd — 2 

(k) 

pairs. Accordingly, here the persistency of entanglement Ve 
reads 

^i''ii*>(^.'^^)i={;i''-iv.,/2£si/2. ■ 

In the perspective of the robustness of the entanglement 
present in the eta pairing states, the direct lattice represen- 
tation seems therefore to be the more favorable, since Ve = 

(k) 

2Ve . More interestingly, the eta pairing states, both in k- 
space and in the direct lattice, display a self-similar behaviour 



under local measures. In particular, once a local measure is 
performed, the output state is characterized by the same struc- 
ture of correlation of the input state, both at the level of bi- 
partite and multipartite entanglement, as at the level of the 
ODLRO present in the state. 

We finally note that recently the possibility of measuring 
correlations between paired momentum modes in a similar 
context have been experimentally achieved, as reported in 

III. BCS STATES 

The direct relation between the entanglement in the mo- 
mentum picture and the ODLRO that has been highlighted for 
the case of the eta pair states suggests to explore whether ana- 
logue relations hold for other states exhibiting superconduct- 
ing correlations. In this last section we thus pass to analyze 
the structure of quantum entanglement in momentum space 
for the first example of state exhibiting superconducting cor- 
relations: the BCS state |2]. We will in particular see how the 
desired direct relation can be established. 
The BCS state was introduced as an Ansatz wave function de- 
fined in momentum representation for models which allow the 
formation of Cooper pairs: 

|BCS) = n(iik + ^k4|aLk|)|vac) =0k|^k) 

k 

(47) 

Here the electron pair is created localized in the pair of modes 
k, — k by the fermionic operators <^k|<^Lk|- The overall state 
is a grand-canonical one since the number of particle is not 
fixed and it is normalized: for each k the coefficients satisfy 
^k + ^k ~ ^- easily see that, at variance with the eta- 

pair states, the structure of correlations is quite simple since 
no multipartite entanglement is present: the state is factor- 
ized in momentum space and the only (quantum) correlations 
exhibited are those between the paired mode k, — k. Indeed, 
the concurrence between the modes k and — k can be evalu- 
ated by treating the |^k) state as effectively qubit states and 
it coincides with the negativity: Ck,-k = 2iik^k = A/k,-k 
1 34] . We now describe how the above entanglement proper- 
ties can be related to the fundamental property of BCS states 
i.e., ODLRO. In order to accomplish this task we can sirnply 
resort to some of the calculations explicitly derived in ll27ll . 
There the two particle density matrix p*^^^ has been obtained 
using the language of Green functions [35] and it has been 
used to study the spin entanglement properties of two elec- 
trons forming a Cooper pair. The explicit expression of p*^^^ 
also allows for the the determination of the ODLRO present 
in the system. We now briefly recall, for the sake of clarity, 
the main ingredients needed for the discussion and we finally 
identify the relevant relations that allow one to link the quan- 
tum properties to the basic superconducting correlations. 
The two-electron space- spin density matrix can be written 

as: 
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(48) 

Here we have that Xi = (xi,^^) and ilj\xi){ip{xi)) is the 
creation (destruction) operator for a particle in position Xi 
with spin Si = ±1/2; the expectation value is taken on the 
ground state at zero temperature. Moreover, the Green func- 
tion G(xiti, ^2^2, x'l^^, ^2^2") = G^(i,2;i',2') IS defined in 
terms of the creation [annihilation] operator 2p\j{xi)[ipH{xi)] 
in the Heisenberg representation as: 

G(l,2;1^20 = (T[^|,(xiti)^|,(x2t2)^i/(44)^i/K^^)]) . 

(49) 

where T is the time-ordering operator and tf is intended as 
a temporal instant following U and infinitesimally close to it. 
Standard arguments [35] show that C(i,2;i',20 can now be can 
be factorized in terms of single particle Green's functions in 
such a way to encompass the presence of pairs of electrons in 
a bound state: 

G(i,2;l',2') = ^^(1,10^(2,20 —G^(l, 20^(2,1') -^(1,2)^(1/ 2') ' 

'(50) 

Here the important part for our discussion is now the anoma- 
lous Green's function: 

i^(\,2) = -i{i^U^ltl)4^H{^2t2)) (51) 

that accounts for the pairing mechanism and is composed by 
the product of a 1/2-spin part Is^s' = 'icfy and a spatial part 
F(xiti, X2ti). The latter is of fundamental importance since 
it is responsible for the ODLRO present in the system; it thus 
is the ingredient in which one can find the connection between 
the entanglement properties of the BCS state in the momen- 
tum representation and its superconducting correlations. In- 
deed, the spatial part of F\ just as much as the spatial part 
of G(^ij), depends only on the differences Axi^2 = xi — X2, 
Si = ti — t[ and it can be evaluated when ^ as: 

iF'''(Axi,2) = ^ ^ UkVk exp [k(xi - X2)] 

k 

= ^IZ-^k,-kexp[k(xi -X2)]. (52) 

k 

The last line comes from the observation that the F^(xi — 
X2 ) can be written as Fourier series whose k components are 
just given by the two mode concurrence previously introduced 
A/k,-k. We finally recall that [27], after taking the continuum 
limit and integrating in the k variable, the limit of the two- 
particles space-spin density matrix when |xi — x-| going to 
infinity defines the ODLRO: 

^^.i,.24;,.^^(xi,2)F*(xi.,20. (53) 

We thus see that, in analogy to what we have described in the 
eta-pair case, in BCS state the superconducting correlations 
that imply Meissner effect and flux quantization are function- 
ally related to the entanglement properties of the state in the 
momentum representation. 



We can finally resume the scenario emerging from the com- 
parison between the structure of correlations of the two kind 
of states considered in this paper. In both cases it possible 
to establish a direct, tough different, link between quantum 
correlations in /c-space and superconducting correlations, see 
(I2TI) and (l52l) . However, while the quantum correlations be- 
tween paired modes play a prominent role in both kind of 
states, the relations found for the eta pairing states between 
the ODLRO and the multipartite correlations suggest that in 
order to establish a superconducting order in lattice models 
that have such states as eigenstates, more resources in terms 
of multipair/multipartite correlations are needed. 



IV. CONCLUSIONS 

In this paper we addressed the task of describing the 
structure of the correlations in the /c-space picture and their 
relations with the superconducting correlations (ODLRO) 
for two important example of states: the eta-pairing and 
the BCS states. These two kind of states share as relevant 
elementary subsystems the paired modes {—kj^kj) and it 
turns out that in both cases the quantum correlations between 
these modes can be directly related with the existing ODLRO. 
In particular, in the eta-paring case the negativity between 
two arbitrary paired modes can be written as the average 
of the pairing correlations that defines ODLRO. As for the 
BCS state, the Green's function formalism developed in ||27|] 
allows to recognize how the ODLRO can be written as a 
Fourier transform of the negativity between the paired modes. 
While the momentum description of the BCS states is fairly 
simple, since the only existing correlations are of bipartite 
nature, the eta-pairing states require a more articulated 
analysis. The latter unveiled a network of correlations which 
is much richer also with respect to the one expressed by the 
same states in the direct lattice picture. Part of the multipartite 
content of the entanglement present in these states is shown 
to have a multi-pair nature and to be directly determined 
by the amount of ODLRO present in the system. The exact 
evaluation of the reduced density matrix of arbitrary large 
blocks of modes enabled the study of the entanglement that 
these subsystems have with the rest of the system (block 
entanglement). In this respect, the peculiar structure of 
correlations in the /c-space, at variance with the direct lattice 
picture, implies that the behaviour of the block entanglement 
in the limit of large numbers of modes heavily depends on the 
peculiar choice of the modes that constitute the subsystem. 
While for blocks composed by Di unpaired modes the block 
entanglement grows linearly with Di, in the case of D2 
paired modes the spectrum of the reduced density matrix is 
given by a binomial distribution, whose variance is given by 
the ODLRO, implying the block entanglement to grow only 
logarithmically with the number of modes D2 . The analysis 
allows one to assess also the intermediate pictures, where an 
arbitrary number D = Di -\- D2 of modes is considered, and 
thus to determine the measure of multipartite entanglement 
firstly devised for qubits by Meyer and Wallach in 1 19]. This 
measure, applied both in direct picture and in the /c-space 
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one, shows that in both cases the multipartite content of the 
entanglement of the eta pairing states is indeed their dominant 
feature and it is thus at the basis of their superconducting 
properties. In view of a possible application to the eta-paring 
states for quantum computational tasks, we finally show how 
the evaluation of the persistency of entanglement points out a 
self- similar structure of the correlations in these states which 
is robust with respect to "local" measurements both in the 
momentum and in direct lattice representation. 
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